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we conduct a general analysis of the whole spectroscopic 
database aimed at identifying and characterizing the 
variety of observed stars, which could be later applied 
to studies of Galactic and extragalactic stellar popu- 
lations. In particular, we define suitable spectral tools 
to determine statistical properties of the stellar sample, 
through their main atmospheric parameters. 



Abstract The GALEX General Data Release 4/5 in- 
cludes 174 spectroscopic tiles, obtained from slitless 
grism observations, for a total of more than 60,000 ul- 
traviolet spectra. We have determined statistical prop- 
erties of the sample of GALEX stars. We have de- 
fined a suitable system of spectroscopic indices, which 
measure the main mid-UV features at the GALEX low 
spectral resolution and we have employed it to deter- 
mine the atmospheric parameters of stars in the range 

4500 < Teff < 9000 K. Our prehminary results in- i xhe GALEX spectroscopic database. 

dicate that the majority of the sample is formed by 
main sequence F- and G-type stars, with metallicity 
[M/H] > -1 dex. 

Keywords ultraviolet: stars; stars: fundamental pa- 
rameters; techniques: spectroscopic 



Introduction 

The ultraviolet (UV) radiation from space is absorbed 
by the Earth atmosphere, thus requiring observations 
from space-borne observatories. At present, the only 
fully UV-dedicated space mission is the Galaxy Evolu- 
tion Explorer (GALEX), which has been observing the 
UV sky since 2003. Up to now, it has provided photom- 
etry of more than 200 million objects (see Bianchi, this 
volume), nevertheless its spectroscopic database is still 
almost completely unexploited. To our knowledge, just 
very few works have been published which made use 
of GALEX spectra (IBello-Ramirez fc Rosa-Gonzalez 



2009HRosa-Gonzalez et al.H2009HBarger fc Cowidboid 
Gowie et al.ll20ld iFeldman et al.ll2010l ). In this work 
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GALEX is a NASA small explorer mission ( Martin et al.^ 
20051 ) whose main goal is to carry out a two-band UV 
photometric all-sky survey. Being equipped with a 
grism, it can also provide slitless spectroscopy of the en- 
tire field-of-view (1.2 deg) in the Far-UV (FUV; 1300- 
1820 A) and Near-UV (NUV; 1820-3004 A), at res- 
olutions of i? = A/AA - 200 and R 100, respec- 
tively. The GALEX General Release 4/^ contains 174 
spectroscopic tiles for a total observing time of about 
2.5 Msec. The database includes the spectra of more 
than 60000 objects. For extracting spectra, we have 
used the automatic pipelin^, which provides the wave- 
length calibrated and flux calibrated monodimcnsional 
spectra of the objects in each tile. 



2 Flux calibration: GALEX vs. lUE 

A reliable flux calibration is undoubtedly a crucial point 
whenever colors are considered or large wavelength in- 
tervals are used to extract stellar atm ospheric param- 
eters. The GALEX flux calibration ( Morrissev et al 



^http://galcx.stsci.edu/GR4/ 

^ |http://www.galex.caltech.edu/DATA/grl_docs/ 
GRl_Pipeline_and_advanced_data_description_v2.htm 
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20051 l2007t) have been obtained from a series of obser- 
va tions of whit e dwar fs (WD) taken from the catalog 
of iBohhn et ahl (|200l[) . To verify the robustness of this 
caUbration, we have compared the WDs in the GALEX 
spectros copic database which are in common with the 
study of iHolberg et ahl (|2003l ). who provide a series of 
high quahty lUE spectra of Galactic white dwarfs. In 
Fig. [U we plot the UV spectral energy distribution 
(SED) of three objects (GD 50, LB 227, and HZ 21) 
that show that the best agreement between lUE and 
GALEX data is achieved in the NUV interval, while 
the GALEX FUV range presents severe flux calibra- 
tion problems, which consist in a overall underestima- 
tion of flux and the appearance of spurious features. 
We therefore limited our present study to NUV wave- 
lengths between 2000 and 2830 A, as, outside this in- 
terval, the GALEX NUV detector sensitivity decreases 
drastically. We would like to mention, however, that 
the "blue" (i.e., short-wavelength) edge of GALEXscn- 
sitivity in the NUV drops at shorter wavelengths than 
lUE, allowing a better spectral morphological analysis 
over the spectral window 2000-2300 A. 



3 The stellar sample. 

Each spectroscopic tile has also been observed in imag- 
ing mode, so that it is possible to assign each spec- 
trum to its corresponding source in the NUV image. 
Wc made use, for each of the 174 spectroscopic tiles, 
of the data present in the photometric catalog file (*- 
xd-mcat.fits) and the spectra provided by the GALEX 
automatic pipeline (*-xg-gsp.fits). Since the goal of this 
work is to assess the potential of GALEX to deter- 
mine the main properties of stars, we extracted point- 
like objects from the entire database by means of the 
following criteria: 

• ellipticity5C0.2 in the NUV image; 

• (SExtractor parameter) class_star^0.8, in the NUV 
image, and 

• a median SNR/px^5 in the NUV spectral interval. 
This criterion by itself decreases the number of spec- 
tra to about 11000. 

The selected analysis sample is formed by a total of 
6037 spectra. 

3.1 Sample properties and non-stellar contamination. 

The adopted selection criteria do not exclude the pres- 
ence in the sample of other kind of sources, besides 
stars. Therefore, in order to have an indication on 
the amount of contamination, we proceeded to cross- 
correlate the position of all sources in our sample with 
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Fig. 1 lUE (Holberg et al. 2002; thick red line) and 
GALEX GR4/5 (thin blue line) spectra of the white dwarfs 
GD50 and HZ21. The GALEX Q.ux error is also shown with 
a dash-dotted line. 

the SIMBAD databascH. We found a total of 2132 ob- 
jects, which is about 35% of our total sample. They 
include 235 entries, which are classified as non-stellar 
(QSOs, galaxies, and so on), and 118 objects which are 
only identified as being sources emitting in UV, X, IR, 
submm, etc. Wc concluded that the expected fraction 
of contamination from non-stellar (or unclassified) ob- 
jects of our sample should be around 16%, which imply 
a total number of stars of about 5000. 

The distribution of the spectral types of the stars 
present in our sample is shown in Fig.O the peak corre- 
sponds to F and G types, with a significant contribution 
of A and K type stars. In an effective temperature scale, 
this translates to roughly the interval 4000-10000 K. 
In this plot we incorporate those stellar objects that 
have a spectral classification in SIMBAD. The fraction 

^ http: / / simbad . u-strasbg.fr / simbad / 
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of white dwarfs with respect of aU other stars can be 
biased from the fact that they were selected for cahbra- 
tion purposes. Examples of high SNR GALEX stellar 
spectra along the spectral type sequence are shown in 
Fig. [3 

The apparent V magnitude distribution (Fig. U) of 
the SIMBAD sample is highly peaked at around 9- 
11 mag. The distribution of stellar objects (thick black 
line) almost coincide with the total distribution up 
to ^15 mag: this implies that they mostly belong to 
the Galactic thi n and thick disk component^ (e-g-, 
Juric et al.ll2008() . Point-like non-stellar objects are the 
majority of dimmer objects in the 15-20 mag interval. 



Table 1 Index passband definition. 



Index name" 


Blue band 


Index band 


Red Band 


gFei2260 


2230- 


-2256 


2259- 


-2285 


2294- 


-2320 


gFeii2332 


2294- 


-2320 


2338- 


-2364 


2421- 


-2447 


gFeii2402 


2294- 


-2320 


2380- 


-2406 


2421- 


-2447 


gBL2538 


2448- 


-2474 


2511- 


-2537 


2644- 


-2670 


gFeii2609 


2448- 


-2474 


2597- 


-2623 


2644- 


-2670 


gBL2730 


2732- 


-2758 






2674- 


-2703 


gMgii2800 


2784- 


-2810 






2674- 


-2703 


g2609/2680 


2597- 


-2623 






2674- 


-2703 


g2000/2180 


2000- 


-2026 






2156- 


-2182 



"The g prefix is added to difi'erentiate these indices 
from those defined in the original Fanelli's (lUE) system. 



*A typical G dwarf of My = -1-5 mag would have a distance of 
158 pc, if my = 11 mag, and 1000 pc, if my = 15 mag. 
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Fig. 2 Spectral type distribution of the subsample of SIM- 
BAD stars. 
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4 A new set of spectroscopic indices for 
GALEX observations. 

The use of spectrophotometric indices to extract se- 
lected information from spectra has been a very com- 
mon pr actice in the optica l; for example, the Lick/IDS 
system (jTrager et al.lll998l and references therein) has 
found a great number of applications in many differ- 
ent research fields. In the UV, however, the use of this 
kind of a nalytical tool has been muc h less widespread. 
Recently, IChavez et all (|2007l . |2009| ) have studied the 
mid-UV spectral morphology o f stars and evolved stel- 
lar populations by means of the lFanelli et al.l (|l992l ) set 
of indices computed in lUE observations. 

For the analysis of our sample, we used a set of 9 
spectrophotometric indices carefully defined to be suit- 
able for the GALEX spectral resolution and to max- 
imize the sensitivity to stellar parameter variations. 
Since the indices mainly depend on the strength of ab- 
sorption features, they are most suitable for late-type 




Fig. 3 Selected GALEX spectra of stars along the spectral 
type sequence. 



Fig. 4 Histogram of the apparent V magnitude of the sam- 
ple of objects with GALEX spectra present in the SIMBAD 
database. The thin red line represents the whole sample, 
while the thick black line indicates the distribution of stel- 
lar objects. 
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Table 2 Linear transformation parameters. 
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Fig. 5 The GALEX spectrum of the star BD+03 40 along 
with the passbands of the set of spectroscopic indices. 



stars, while they decrease their effectiveness for hot O- 
B stars. 

The index ID and passbands arc reported in Table [T] 
and shown in Fig. [5l For 5 i ndices we have ado pted 
the definition given by Eq. 5 in lFanelli et al.l (|l992l ). In 
this definition the flux in the index band is compared 
to a pseudo-continuum, obtained by the interpolation 
of the fluxes in the two side bands; the other 4 indices 
are narrow-band colors, which measure the difference 
in magnitude between the blue and red band. Note 
that, due to the red edge cut-off sensitivity, it was not 
possible to define a suitable red band for the indices 
measuring the strength of the Mg ii doublet at 2800 A 
and the blend around 2730 A: we have therefore defined 
two colors, where the blue and red bands are inverted; 
in this way, the indices have increasing values with in- 
creasing line absorption. The first 8 indices of the table 
measure the intensity of the same spectral features of 
the Fanelli et al. system, apart from gBL2730, which 
merges together the two lUE absorption blends at 2720 
and 2740 A due to the lower GALEX spectral reso- 
lution. The brand new index that we have defined, 
the color g2000/2180, takes advantage of the higher 
sensitivity of GALEX at A < 2300 A with respect to 
lUE. This property also allows more precise measure- 
ments of the gFeil2 332 and gFeii2402 indices, which 
Chavez et al. ( 2009f ) point out as being good tracers of 
relatively young and metal rich stellar populations. 

The new indices were measured in three different 
spectral samples: 



• the 6037 selected GALEX point-like sources; 

• the ATLAS library of synthetic stellar SEDs of 
Castelli &: Kuruc"3 ( 2003H ^. after having homogenized 



Index 


a 


b 


<j 


used 


excl. 


gFei2260 


0.016 


0.597 


0.022 


71 


23 


gFeii2332 


0.029 


1.369 


0.128 


79 


15 


gFeii2402 


0.062 


1.699 


0.177 


82 


13 


gBL2538 


0.025 


1.419 


0.095 


61 


34 


gFeii2609 


-0.008 


1.285 


0.097 


67 


28 


gBL2730 


0.023 


1.430 


0.127 


79 


15 


gMgii2800 


0.223 


1.133 


0.274 


84 


10 


g2609/2680 


0.019 


1.151 


0.088 


72 


23 


g2000/2180 


0.125 


1.365 


0.426 


80 


13 



them to the same spectral resolution and wavelength 
sampling of the GALEX data. Since the flux is 
sampled at AA = 10 A, this library is suitable for 
analysing GALEX spectra; 



a cho sen sample of lUE stars, from the IWu et al 



(|1983[ ) catalog, whose atmospheric parameters are 
known. In this case, we have also properly degraded 
the original lUE spectra to the GALEX resolution. 
This set of objects was used in the calibration process 
described in the following Sec. 14.11 



4.1 Calibration of the theoretical indices. 



Since the ATLAS SEDs do not perfectly rep r oduce 



the spectra of real sta.rs (IBertone et al.l l2004 12008 



Rodriguez-Merino et al.l l2005l ) . the theoretical indices 



must be converted to the observed system, before ap- 
plying them to the analysis of the GALEX sarn p le. 



performed this calibration, as in IChavez et al.l (|2007() , 
by means of a sample of stars, which were observed 
by GALEX and whose complete set of atmospheric 
paramet ers (Tf.ff,log (7,[M/Hl)_ is given in the PASTEL 



catalog ( Soubiran et al. 2010t ). Since the total num- 



http : / / w wwuser . oat. ts. astro. it/castcUi/ 



ber of these stars is not large enough to provide a 
goo d calibration, we added a sample of lUE star from 
the IWu et al.l (|l983l ) catalog. The theoretical spectra, 
with parameters corresponding to each of the calibra- 
tion s tars, were obtained by a tri-linear interpolation 
of the ICastelli fc Kuruczl ([2003!) grid. Figure |6| shows 
the comparison between the observed and the theo- 
retical indices for three indices, which are representa- 
tive of the different degrees of agreement, from one of 
the best cases (g2609/2680) to the worst (gMgii2800). 
The trasformation equations are obtained by an it- 
erative process of linear interpolation (indexthoor = 
a + b X indexobs), excluding at each step the outliers 
which stand more than 2a away from the fitted value. 
In Tab.|2]we report, along with the a, 6 parameters, the 
final rms and the number of stars used and excluded 
in the fitting process. There is a clear tendency of the 
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g2609/2680 gFell2402 gMgll2800 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 -0.6-0.4-0.2 0.0 0.2 0.4 0.6 0.0 0.2 0.4 0.6 0.8 1.0 1.2 

observed (GALEX + lUE) observed (GALEX + lUE) observed (GALEX + lUE) 



Fig. 6 Theoretical vs. observed indices for tlie GALEX (blue circles) and lUE (red triangles) samples. The solid black 
lines indicates the best linear fit, which did not take into account the outliers (open symbols); the parallel dotted lines lie 
fcr apart. The green dashed line shows the one-to-one correlation. 



theoretical value to be higher than the observed data, 
the only exception being the gFei2260 index. This dis- 
agreement makes the tranformation process compulsory 
if the goal is to apply this set of theoretical indices to the 
determination of the atmospheric parameters of stars. 

4.2 Theoretical indices vs. stellar atmospheric 
parameters. 

Within the theoretical framework, we explored the sen- 
sitivity of each index with respect to the fundamen- 
tal stellar atmospheric parameters: effective temper- 
ature (Toff), surface gravity (logg), global metallicity 
([M/H]), and the [a/Fe] enhancement. The results are 
displayed in Fig. [7l The indi ces that measure co mmon 
spectral features with those of lFanelli et al. I (ll992l) show 
similar trends, with respect to Toff, logg, and [M/H], as 
those calculated with the Fanelli's passbands definition 
( Rodriguez- Merinol[200l IChavez et all l2007t ). Apart 
from the two bluest indices, the values are larger for 
lower gravities. The slope g2000/2180 is an interest- 
ing case, since it is almost independent of gravity over 
the whole temperature range. As far as [M/H] is con- 
cerned, all indices show a strong dependence and a sim- 
ilar pattern: the temperature at the peak increases with 
in creasing metal abundance. 

Chavez et ail |2009f ) found a strong effect of the a- 



enhancement in the SEDs of stars and simple stellar 
populations all over the mid-UV interval. In what fol- 
lows, we explore its impact on t he spectroscopic indices. 
In the Castelli fc Kuruc"3 ( 2003 ) library, the a-enhanced 
models and fluxes were computed by adding an extra 
abundance of -1-0.4 dex to each a element, which causes 
an increase of about -1-0.3 dex to the global metallicity 
Z with respect to the nominal one. Therefore, to com- 
pare the results at equal Z, we scaled down all abun- 
dances of the [a/Fe] =-1-0.4 models to reach the same 



global metallicity as the solar-scaled models and we 
computed the correspondent SEDs by linear interpo- 
lation in the grid of theoretical fluxes. The abundance 
rescaling produces a significant reduction of the non-a 
elements, in particular of the Fe, which is reflected in 
a decrement of almost all the GALEX indices, whose 
intensity is strongly dependent from the absorption of 
iron atoms, which is one of the main sources of opacity 
in the mid-UV. The indices gBL2538 and gMgii2800 
are the only ones which present non negligible posi- 
tive differences in a-enhanced spectra: this is due, in 
both cases, to the presence of magnesium and, subor- 
dinately, titanium absorption lines inside the central 
index band. In Fig. [71 we also tra ce a comparison of 
the in dices computed for the original lCastelli fc Kurucz 
( 20031) [a/Fe]=-h0.4 SEDs and the correspondent solar- 
scaled ones: in this case, besides the effect of the en- 
hancement of the a elements, the difference takes into 
account the global increasing in metallicity; therefore, 
the difference has higher values (i.e., less negative) at 
higher temperatures, where the index strength increases 
with increasing global metallicity, and, vicevcrsa, a- 
enhanced indices are lower were the index dependence 
to metallicity is reversed. 

We conducted a complementary analysis assessing 
the effect of interstellar extinction on the indices. For 
this test we considered E{B — V)=OA. which is a typical 
value for the objects of our sample. The results indicate 
that index variations due to extinction arc negligible. 



5 Atmospheric parameters of the GALEX 
objects. 

We have used the set of indices defined here to deter- 
mine the three main stellar atmospheric parameters for 
all the objects of the selected sample. For non-stellar 
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Fig. 7 Left column: Theoretical indices, at fixed [M/H]=0.0, as a function of effective temperature for different values of 
log 3= 2 (red curve), 3 (yellow), 4 (green), and 5 dex (blue). Central column: Theoretical indices, at fixed logg=4 dex, as a 
function of effective temperature for different values of solar-scaled metallicity [M/H]= -2 (red curve), -1 (yellow), (green), 
and +0.5 (blue). Right column: Effect of the a-enhancement on the theoretical indices for the same metallicity values 
shown in the central column: the thick curves show the [a/Fe]=+0.4 minus the solar-scaled indices once the abundances 
are re-scaled to have the same Z as the solar-scaled models, while the thin dashed lines show the same difference obtained 
by adding 0.4 dex to the number density of a elements. 
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Fig. 8 xit contour plots in tlie logp vs. Tcs plane (upper 
panels) and [M/H] vs. T^b plane (lower panels) for 3 GALEX 
objects. Tlie filled circles mark the loci of the minimum x^, 
while the grey surfaces bounded by the 1-, 2-, and 3-(t er- 
ror contours. In the panels are also indicated the GALEX 
identification code, the best values of the atmospheric pa- 
rameters, and of the minimum x^- 



objects contaminating the sample the results are, of 
course, meaningless. In this first work on G^LiJX spec- 
tra, we restricted the examination of the results to a 
statistical analysis. 

We obtained effective temperature, surface gravity, 
and global metallicity of each object by minimizing the 
chi-square distances between the observed set of indices 
and the theoretical ones. The reduced chi-square is 



defined as: 



2 1 V (-^thcor.ij 



V 



a, 



(1) 



obs,i 



where i runs over the spectroscopic indices /, while j 
indicates the theoretical model and v is the number of 
indices which could be computed from the spectrum. 
The error a for each observed index is determined via a 
Montecarlo computation; a random error taken from a 
Gaussian distribution, with standard deviation equal to 
the observed flux error, is added to the original GALEX 
spectrum and the indices are computed on this new 
spectrum. This is repeated 200 times and the error of 
each index is given by the standard deviation of the 
index distribution. 

We searched for the combination of Toff, log 5, and 
[M/H] that presents the minimum x^, so that we could 
extract the planes at fixed metallicity and at fixed grav- 
ity that included this point. We then refined the de- 
termination of the minimum on the two surfaces after 
smoothing them with a quintic surface, by means of 
the tri_surfliyL function. The loci of the minima allow 
the determination of the best values of Tog, logg, and 
[M/H] of each object. In Fig. HI we present three exam- 
ples of the contour plots of the xt in the \ogg vs. Toff 
and [M/H] vs. T^s planes: they show one case where the 
procedure provides a precise determination of the stel- 
lar parameters (upper panel), a second case where the 
roff-[M/H] degeneracy generates two l-cr regions (mid- 
dle panel), and finally an example of a failure of the 
method to determine the parameters with an adequate 
precision (lower panel). 

The distributions of the results are shown in Fig. HI 
The histogram of the effective temperature peaks at 
TcS ~ 6000 K, typical of carly-G or late-F spectral 
types, and show very large numbers up to Tcff ^ 7000 K 
(the F-type interval). A significant fraction of objects 
have a temperature of later G- and K- (approximately 
4000-5800 K) and A-type (-7000-10000 K) stars, while 
outside this interval the numbers are very small. This 
distribution well agrees with the stellar type histogram 
of the SIMBAD stellar sub-sample (Fig.0). 

The surface gravity distribution indicates a prepon- 
derance of main sequence objects (logg > 4 dex) with 
respect to more evolved stars, while the histogram of 
the global metallicity shows a majority of stars with 
[M/H] > — 1 dex. This is in agreement with the prop- 
erties of the Galacti c stellar populatio ns of the thin 
and thick disks (e.g., iBuser et al.lll999i and references 
therein). 

Since the theoretical stellar library includes only two 
values of a-enhancement, we excluded, for now, the 
determination of this further parameter. We can say. 
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Fig. 9 Parameter distributions of tiie GALEX sample. 
The thick black lines refer to the whole sample, while the 
thin red lines show the histograms of the results which have 
a significance level greater than 1%. 



nevertheless, that, according to the behaviour of the 
indices with [a/Fe] shown in Fig. [7J we expect to have 
slightly overestimated the effective temperature and un- 
derestimated the global metallicity if an enhancement 
is present in some of the stars. 
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